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ABSTRACT 

Peculiar velocities throughout the region of the Local Supercluster are reconstructed by two different 
orbit-retracing methods. The requirement of the optimal correlation between the radial components 
of reconstructed velocities and the observed peculiar velocities derived from our extensive new catalog 
of distances puts stringent constraints on the values of the cosmological parameters. Our constraints 
intersect those from studies of microwave background fluctuations and statistical properties of galaxy 
clustering: the ensemble of constraints are consistent with fi m = 0.22 ±0.02. While motions through- 
out the Local Supercluster provide a measure of the mean ratio of mass to light, there can be large 
local fluctuations. Our reconstruction of the infall velocities in the immediate vicinity of the Virgo 
Cluster shows that there is a mass-to-light anomaly of a factor of 3 to 6 between groups in the general 
field environment and the heavily populated Virgo Cluster. 

Subject headings: dark matter — cosmological parameters — methods:analytical and numerical 



1. INTRODUCTION 

Cosmological parameters have been significantly con- 
strained by the combined analyses of the data mainly 
from Cosmic Microwave Background (CMB) radiation, 
large-scale galaxy distribution, supernova observation s 
and cosmic shear field [see e.g. l|Tegmark et alJ l20041]. 
Here, we demonstrate that analysis of the peculiar ve- 
locity field can tighten the constraints on cosmological 
models because error ellipses are strongly skewed to those 
of WMAP and SDSS in the parameter space of density 
and age. Two complementary theoretical methods are 
used in this paper to generate orbits that give veloci- 
ties to be compared with observations. Peculiar veloc- 
ities are defined by a large number of quality distance 
measures. Good distances are available for 1400 galaxies 
within 3,000 km s _1 . 

The two variational methods used here are based on 
the principle that galaxies as mass-tracers follow orbits 
which are the stationary points of the Eulcr-Lagrange 
action. Our fir st approach, the Least Action (LA) method 
l|Peebleslll989D , describes orbits in the nonlinear collapse 
regime, though recovered trajectories are non-unique due 
to lack of knowledge of the initial positions of galaxies 
and because of orbit-crossing (i.e. multisteaming) . On 
large scales, a unique reconstruction is possible since the 
displacement of dark matter fluid derives from a potential 
and multistreaming is less severe. Th ese properties are a t 
the basis of our second method, MAK ( Fri sch et all 2002') . 
The methods are applied at the Local Supercluster scale 
for the determination of the cosmological mean density 
and at the cluster scale for the specification of its local 
variations. 

Allowing for statistical and systematic errors and con- 
sidering all the methods in tandem, we show that the 
MAK, LA, WMAP and SDSS results give preference to a 
small region in parameter space favouring a low value of 
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f2 m = 0.22 ± 0.02 for the density parameter. 

On cluster scales, reconstruction of the amplitude of 
the first-approach infall velocities provides a determina- 
tion of the mass of the clusters at radii beyond the viri- 
alized region; radii greater than those explored by X-ray 
or lensing studies. We show that the high amplitude 
observed in the infall region requires an assignment of 
M/L = 5OOM /L to the Virgo Cluster (where the lu- 
minosities are in the B band). 

2. CATALOG OF NEARBY GALAXIES AND DISTANCES 

The catalog of galaxies, is a 40% aug mentation of 
the Nearby Galaxies Catalog l)Tullvlll987|) . now includ- 
ing 3300 galaxies within 3,000 km s _1 . This depth is 
more than twice the distance of the dominant compo- 
nent, the Virgo Cluster, and the completion to this depth 
in the current catalog comp ares favorably with other 
all-sky surveys [e.g., 2MASS ij.Tarretd 120041) ]. The cat- 
alog has the following properties, (i) All entries are 
given a detailed group and filament assignment (or non- 
assignment if isolated), (ii) The catalog is supplemented 
by an all-sky complet e sample of X-ray selected clusters 
l|Kocevski et al.ll2004|) in the shell 3,000-8,000 km s" 1 to 
provide a description of potential influences on very large 
scales, (iii) 'Fake' galaxies are added at low Galactic lati- 
tudes to avoid an underdensity in the zone of obscuration 
due to lost information, (iv) Correction is made for the 
loss of light with distance caused by an apparent magni- 
tude cutoff in the construction of the catalog. 

The second observational component is an extended 
catalog of galaxy distances. Information from four 
techniques has been i ntegrated: the Cepheid variable 
( Freecjman et^ al.||2 1[). Tip of the R ed Giant Branch 
( Karachentsevl 12003 iLee et alJ 119931) Surface Bright- 
ness Fluctuation ijTonrv fc Schneiderlll988t iTonrv et all 
20011). and Lum i nosity -Linewidth l|Tullv fc Fisherlll977t 
Tull v & Pearcel [2000) methods. In all, there are 
over 1400 galaxies with distance measures within the 
3,000 km s _1 volume; over 400 of these are derived by 
at least one of the first three 'high quality' techniques. 
In the present study, distances are averaged over groups 
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because orbits cannot meaningfully be recovered on sub- 
group scales. The present catalog is assembled into 1234 
groups (including groups of one) of which 633 have mea- 
sured distances. These observational components will be 
described in detail in a later publication. 

3. TECHNIQUES: THE LA AND MAK METHODS 

Peebles pioneered the Least Action (LA) method 
<|PeehlesHl98fl 179941 ITflfll that requires orbits to sat- 
isfy the stationary point of the Euler-Lagrange action, 
the integral over time of the Lagrangian. Hence, least 
action searches for the minimum of 
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S LA =/ dt -x _ + Gph0 2 m . x 2 (1) 
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where summation over repeated indices and j ^ i is im- 
plied, to denotes the present time, the path of the ith 
particle with mass rrii is Xj(i), pb is the mean mass den- 
sity, and the present value of the expansion parameter 
a(t) is ao = a(to) = 1. 

Individual orbits are constrained by the mixed bound- 
ary conditions that peculiar velocities were initially negli- 
gible (they subsequently grew from gravitational pertur- 
bations) and the known elements of position and velocity 
today. The other piece of information known for roughly 
half the elements is the aforementioned distances. This 
information is used to discriminate between models. A 
given model involves a specification of cosmological pa- 
rameters and the assignment of mass to each of the ele- 
ments. Orbits are found within the context of a specified 
model that are consistent with the boundary condition 
constraints on angular positions and radial velocities. A 
component of the end point of an orbit is its distance 
which can be compared with the observed distance. Dif- 
ferent models result in different distances. The quality of 
a model can be evaluated by a statistical measure of the 
differences between observed and model distance mod- 
uli, ^obs and \i mo d respectively. A \ 2 parameter can be 
calculated for each element, i, with a measured distance, 
Xi = (fJ-mod.i - fJ-obs,i) 2 /£i, where ti is the uncertainty 
assigned to [i b s ,i- A comparison of the relative distri- 
bution of xf values provides a sensitive discriminant be- 
tween good and bad models ijShava et al.lll995t iPhelpgl 
2002). 

Galaxy flows a re not sensitive to any uniform repu l- 
sive dark energy ijLahav et al. 1 11991 : Sh ava et al.l ll995h 
Measurements of the lowest frequency peak in the mi- 
crowave background fluctuation spectrum strongly indi- 
cate, and we accept, that the Universe has a flat topol- 
ogy, whence Q\ + £l m = 1 where Q\ characterizes the 
mean density of dark energy and Q m characterizes the 
mean density of matter with respect to the critical den- 
sity for a flat Universe in matter alone. Specification of 
Q m and the Hubble Constant, H , which describes the 
mean expansion rate of the Universe, in the context of a 
flat Universe uniquely specifies the age of the Universe, 
io(^mj Ho). 

The Least Action orbits are a reflection of the inter- 
actions of masses over the age of the Universe. As an 
initial approximation, it is assumed that the relationship 
between mass and light is constant so the mass Mi as- 
signed to an element in the catalog with luminosity i,; is 
Mi = (M/L)Li. Consequently, a cosmological model to 
be explored is specified by two fundamental parameters: 



to(£l m , Ho) and M/L. For a given choice of to and M/L, 
orbits are constructed, model distances at the end points 
are determined and compared with observed distances, 
and the match provides a \ 2 measure that can be com- 
pared with alternative models. Discrimination between 
results over the domain of (to, M/ L) choices based on the 
X 2 measure gives specification of a best model. 

The alternate technique that will be implemented 
here is the Monge-Ampere-Kantorovich (MAK) re- 
construction jFrisch et alJ 120021: iMohavaee et all 120031: 
iBrenier et al.ll2003[) . This method provides a recipe for 
the orbits that is unique to the degree that orbits can 
be described as following straight lines under suitable 
coordinate transformations. Orbits in the MAK recon- 
struction are minima of the action which assume a La- 
grangian mapping qnx that can be described in terms 
of a potential as x = y q $(q). Here, x is the current 
Eulerian position and q is the initial Lagrangian posi- 
tion. The potential $(q) is assumed to be convex; ie, 
orbit-crossing is excluded. With these conditions, orbits 
are reconstructed uniquely by an assignment algorithm 
that finds the minimum in the ensemble of orbit dis- 
tances summed in quadrature. Hence, MAK search for 
the minimum of the discretized action 
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Initial positions , are assigned to final positions x^ . For 
each final position i, there are N possible initial positions 
and only the one that minimizes the sum in the expres- 
sion (2) is allowed. A brute-force assignment algorithm 
would have a complexity of N\, but can be reduce to a 
polynomial complexity (specifically, so far to TV 2 25 ). It 
is insightful to note that the least action variation Q 
reduces to MAK optimization © for inertial trajectories 
(i.e. where particles move with their initial velocities). 
The distances, d, permit an extraction of peculiar ve- 



locities V m 



v a! 



dHo where V qsr is the observed 



velocity of an object in the galactic standard of rest. 
The relationship between redshift and real space dis- 
placements of the elements is estimated using the ap- 
proxi mation that galax ies trace mass on inertial trajec- 
tories l)Zerdovichl ll970') : v = /H(x — q) where v is the 
peculiar velocity vector and / ~ ^m 6 - Consequently, a 
specific model defines positions that can be tested against 
observed positions in the same way that is done in the 
Least Action analysis. 

An advantage of Least Action is that it potentially 
can be used in the highly non-linear regime where orbits 
are strongly curved. An advantage of MAK is that it 
recovers orbits uniquely and quickly, hence can be used 
with the very large catalogs that are becoming available. 
By using the two methods on the identical catalogs, with 
the identical distance constraints, it is affirmed that the 
two methods recover similar results in the regimes where 
they are both applicable. 

Both methods have been tested on N-body simula- 
tions and found to result in modest underestimations 
of the known densities. Underestimates are anticipated 
because, at one extreme, mass that is not strongly 
clustered on the scale of the survey has little dynam- 
ical consequence and, at the other extreme, complex 
shell-crossing orbits cannot be modeled. Analysis of 
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Fig. 1. — Reconstruction methods constrained by precision dis- 
tance measures give estimates of the mean mass density of the 
Universe. The contours are \ 2 a t 1 an d 2<r for models with dif- 
ferent choices of age and mass density Qm, with shaded contours 
for MAK and line contours for LA reconstructions. The MAK 
and Least Action contours are adjusted upward in density by 20% 
and 30% respectively to compensate for systematic underestima- 
tions of density. The confluen ce of the cons traints on density and 
age paramet ers from WMAP ^ pcrgcl 2003) of Q m h 2 = 0.134 and 
from SDSS ITegmark et al.120041) of h = 0.21 is lightly shaded. 
The darkest shaded region lies within the 2a limits of all four ex- 
periments. The right panel shows the density determinations for 
LA, MAK, SDSS, and WMAP for H = 80 km s" 1 Mpc' 1 , the 
value of the Hubble Constant consistent with our distance scale 
zero point. The open and filled symbols for LA and MAK respec- 
tively give the results before and after compensation for the system- 
atic under-determination of mass density. At H = 80, t$ bss = 



12.2 Gyr, 



analyses are consistent with £7 r , 
value of Hq . 



13.0 Gyr, and t 



13.5 Gyr. The 4 
0.22 ± 0.02 with the constrained 



MAK reconstructions in 12 cosmological simulations with 
different initial fluctu ation characteristics discussed by 
iMohavaee et alJ (|2005f) resulted in recovery of 80% ±20% 
(la) of the known model density. In the case of Least 
Action, it has been appreci ated that mass can be under- 
estimated on small s cales (jBranchini fc Carlberd IT"994t 
iBranchini et al.lf2002fl . A quantitative measure of the ef- 
fect is found using our Least Action reconstruction al- 
gorithms with an N-body simulation with Q\ = 0.7, 
fi TO = 0.3, evaluated by placing the observer at mul- 
tiple locations to test for cosmic variance. Recovered 
densities were 70% ± 20% of the model density (Phelps, 
Desjacques, and Nusser; ongoing work). The estimates 
that we quote for fi m include adjustments for the sys- 
tematics of 20% and 30% for MAK and Least Action, 
respectively. 

4. RESULTS: I. DETERMINATION OF Q M 

The results of the two methods are summarized in 
Fig. 1 where contours of the % 2 parameter are plotted 
in the domain {to,il m }. In the case of the Least Ac- 
tion analysis, mass density Q m is derived from M/L 
values by accepting the mean B band l uminosity den- 
sity of the Universe determined from SDSS lBlanton et aP 
(2003) with reddening corrections. We use fl m = 
(M/L)/(1540/i), a transform with - 20% uncertainty. 

The error ellipses from the MAK and LA studies are 
elongated. With a linear analysis the elliptical \ 2 troughs 
would open to infinity, but non-linear effects create a 
specific minimum along each \ 2 trough. At a fixed age, 
to, there are relatively tight constraints on Q m require- 



ments. If shorter ages are entertained, then higher den- 
sities are required to arrive at the observed dynamical 
state in the specified time. The contrary dependence 
between density and time are seen in the results from 
the microwave backg round WMAP (iSpergell 120031) and 
galaxy redshift SDSS l)Tegmark et al.ll2004j) experiments 
that are also superimposed on Fig. 1. Only a small do- 
main around to — 13.5 Gyr lies within or near the 2a 
contours of all the methodologies. With the constraint 
h = 0.8 consistent with the zero point of the distance 
estimates, there is good agreement between the WMAP, 
SDSS, MAK, and LA measures of the density parameter: 
TO = 0.22 ± 0.02. The low value of O m obtained here 
is consistent with that found from divers recent studies 
(iBahcall et al. l2003HOstrikerl200llvan den Bosch et alJ 
l2003UCole et al.ll2005|) . ~ 

5. RESULTS: II. THE MASS OF VIRGO CLUSTER 

The velocity field analyses discussed above make the 
simplistic assumption of constant M/L. There has been 
evidence in the literature, t hough, for large variations 



in M/L with environment llMarinoni fc Hudso 



van den Bosch et al 



2003: iFke et nl.ll20l. lt iHTull 
2005; IParker et a,Ul2005jL I 
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2002; 



2005; 



van den Bosch et a l. 2005; Park er et alJ l2005L Here we 
note that galaxy flows near the Virgo Cluster contra- 
dict the hypothesis of constant M/L. The high veloci- 
ties of infall suggest a much larger M/L for the cluster 
than the value consistent with a good fit over the full 
region within 3,000 km s _1 . The critical galaxies for 
this discussion lie outside the cluster on the plane of the 
sky so only a modest component of the infall motion is 
projected into the line of sight but there is no confu- 
sion with cluster membership. Galaxies infalling from 
the foreground of the cluster are redshifted with respect 
to the cluster and galaxies infalling from the background 
are blueshifted. If one could imagine many test parti- 
cles distributed along a line of sight, they would create a 
wave in velocity as a function of d istance with the same 
velocity arising at three distances ijTonrv &: Davisiri981|) 
- two located within the infall region and one located 
at the cosmic expansion position. The amplitude of the 
peaks of the 'triple-value' waves depend on the mass, 
M, interior to the position associated with the peaks, r. 
The envelope of observed infall vel ocities with r provides 
a description <|Tullv fc Shaval 19841) of the run of M(< r). 
The outermost caustic of the Virgo Cluster proper (the 
radius of second turnaround) is at ~ 2 Mpc from the 
center of the cluster. Galaxies on first infall acquire very 
high radial velocities at 2-4 Mpc from the Virgo core. 

An example is given in Fig. 2. The run of velocities 
with distance are what is expected along the specific line 
of sight of the 11-3 group in the Nearby Galaxies Catalog 
with three different assumptions regarding the mass of 
the Virgo Cluster. The observed distance and velocity 
of the 11-3 group is indicated in the figure. It is clear 
that the 11-3 group is at one of the two infall locations 
(or the location where the two infall positions become 
degenerate at the tip of the triple-value curve) rather 
than at the cosmic expansion position around 30 Mpc. 

The total mass within 3,000 km s _1 is the same in all 
3 models to ensure that the global % 2 is near the mini- 
mum of the trough in Fig. 1, however the distribution of 
mass between the elliptical rich and spiral rich groups is 
modified between models. Here, to = 13 Gyr. In the case 
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Fig. 2. — The distortion of the velocity-distance relation by 
the Virgo Cluster in the line of sight toward the 11-3 group. The 
observed distance and velocity of the 11—3 group is indicated by 
the point with error bars in distance (errors in velocity are smaller 
than the symbol) . The observed distance and velocity of the Virgo 
Cluster is indicated (centered at 6.2 degrees from the line of sight 
of the 11—3 group). The background Virgo W Cluster causes the 
second wiggle along this line of sight. The behavior of infall into 
the Virgo Cluster is traced by the locus of positions and velocities 
of test particles of negligible mass that were laid down in the line of 
sight of the 11—3 group and subjected to Least Action modeling. It 
is concluded from the information seen on this plot and equivalent 
plots of the lines of sight to other groups close to the Virgo Cluster 
that M/Lvirgo ~ 5OOM Q /L0 and the mean M/L for the field 
outside the cluster is approximately 1/3 of this value. 



of the curve with the smallest swing, all components are 
given the same M/L = 200M©/ Lq. The intermediate 
curve is generated with the elliptical-dominated groups 
(including the Virgo Cluster) given M/L E = 500 and 
the spiral-dominated groups given M/Ls = 174. The 
curve with the largest swing reflects M/Le = 800 and 
M/Ls = 148. In the model illustrated by the middle 
curve in Fig. 2, the cluster is given the necessary and 
just sufficient mass of 9 x 10 14 M Q . 

The 11-3 group and a couple of others provide the 
greatest demands on the mass of the cluster and are in- 
ferred to lie near a tip along a triple-value curve. An 
adequate description of the amplitude of their infall ve- 
locities requires an assignment of M/L = 500Mo/L© to 
the Virgo Cluster if to — 13 Gyr. As with the overall Su- 
percluster modeling, less cluster mass is required if given 



more time llTullv fc Shavalll98l : e.g., 8 x 10 14 M Q suf- 
fices if to = 14 Gyr. In both these 13 and 14 Gyr cases, 
as the mass assigned to the cluster is augmented, the 
models require a reduction of mass assigned to the field 
to remain at the minimum of the \ 2 trough of Fig. 1. 
The ratio of cluster, C, to field, F, mass to light val- 
ues is (M/L)c/(M/L)p — 3 with the mass assignments 
required to explain the infall motions. This factor 3 prob- 
ably underestimates the M/L difference between Virgo 
and bound groups in the field because some of the mass 
contributing to the field ratio lies outside the groups . 
Studies of the dynamics of nearby groups l|Tullvl 12005^ 
suggest that bound groups of spiral galaxies in the field 
have M/L ~ 90. The mass to light ratio of the Virgo 
Cluster is 5-6 times larger. 



6. CONCLUSIONS 

Our first result on the mean density of the Universe 
complements the WMAP and SDSS measurements of the 
density parameter because the error ellipses in the do- 
main Q m ,Ho are steeply inclined to each other. A nar- 
row range of parameter space around tg = 13.5 Gyr, 
f2 m = 0.22 ± 0.02 is permitted by the combined experi- 
ments. 

Our second result concerns the local fluctuations 
around this mean density. Although most of the blue 
light in the local Universe is associated with the spiral 
field, there is close to parity in the partition of mass be- 
tween the spiral field and the Virgo Cluster. The M/L 
value associated with the E/S0 dominated Virgo Clus- 
ter is significantly higher than that associated with the 
spiral dominated field. The complex problem of M/L 
variations with environment shall be discussed in a fu- 
ture more-detailed paper. 
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